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In situ spacecraft data on the solar wind show events identified as magnetic reconnection with 
outflows and apparent “X-lines” 10 3-4 times ion scales. To understand the role of turbulence at 
these scales, we make a case study of an inertial-range reconnection event in a magnetohydrodynamic 
(MHD) simulation. We observe stochastic wandering of field-lines in space, breakdown of standard 
magnetic flux-freezing due to Richardson dispersion, and a broadened reconnection zone containing 
many current sheets. The coarse-grain magnetic geometry is like large-scale reconnection in the 
solar wind, however, with a hyperbolic flux-tube or “X-line” extending over integral length-scales. 


Magnetic reconnection is widely theorized to be the 
source of explosive energy release in diverse astrophysi- 
cal systems, including solar flares and coronal mass ejec¬ 
tions [1], gamma-ray bursts [2], and magnetar giant flares 
[3]. Because of the large length-scales involved and con¬ 
sequent high Reynolds numbers, many of these phenom¬ 
ena are expected to occur in a turbulent environment, 
which profoundly alters the nature of reconnection [4-6] . 
In the solar wind near 1 AU, which is the best-studied 
turbulent plasma in nature, quasi-stationary reconnec¬ 
tion has been observed for magnetic structures at a wide 
range of scales, from micro-reconnection events at the 
scale of the ion gyroradius (^100 km), up to integral 
length scales 10 4-5 km), and even to larger scales 
[7] . Yet numerical studies of reconnection in magnetohy¬ 
drodynamic (MHD) turbulence simulations have focused 
almost exclusively on small-scale reconnection at the re¬ 
sistive scale [8-10]. Our objective in this Letter is to 
identify an inertial-range reconnection event in an MHD 
turbulence simulation and to determine its characteristic 
signatures, for comparison with observations in the solar 
wind and other turbulent astrophysical environments. 

To search for reconnection at inertial-range scales we 
adapt standard observational criteria employed for the 
solar wind. In pioneering studies, Gosling [] ] has looked 
for simultaneous large increments of magnetic field 5B(r) 
and of velocity field 5u(r) across space-separations r near 
the proton gyroradius p p , which approximate MHD ro¬ 
tational discontinuities. Candidate reconnection events 
are then identified as pairs of such near-discontinuities, 
with 5B(r) aligned for the two members of the pair and 
5u(r) anti-aligned. Gosling’s selected events generally 
have the appearance of two back-to-back shocks, or a 
“bifurcated current sheet”. We modify this criterion to 
allow for more gradual field-reversals, by choosing instead 
r = L/10, with L the outer (integral) length of the tur¬ 
bulent inertial range and by considering pairs separated 


by distances up to L/ 2. 

We apply the above criterion to two datasets. The 
first is from a numerical simulation of incompressible, re¬ 
sistive MHD in a [—7r,7r] 3 periodic cube, in a state of 
stationary turbulence driven by a large-scale body-force. 
The simulation has about a decade of power-law inertial- 
range and the full output for a large-scale eddy turnover 
time is archived in an online, web-accessible database 
[11]. The second dataset consists of Wind spacecraft ob¬ 
servations of the solar wind magnetic field B, velocity 
u, and proton number density n p . The results presented 
here are from a week-long fast stream in days 14-21 of 
2008 (cf. [12]). The average solar wind conditions were 
u = 638 km/s, B = 4.3 nT, n p = 2.3 cm -3 , Alfven speed 
Va = 62 km/s, proton beta /3 P = 1.1, and proton gyrora¬ 
dius p p = 154 km. The temporal data-stream from the 
spacecraft is converted to an equivalent space series using 
Taylor’s hypothesis, x = ut [13]. Simulated “spacecraft 
observations” from the MHD database are taken along 
192 linear cuts, with 4 3 cuts through each face of the sim¬ 
ulation cube in the three coordinate directions. We find 
a good correspondence for statistics of 5B(r) and 5u(r) 
in the two datasets, with the grid spacing dx = 27r/1024 
of the simulation related to 18 s of the Wind time-series 
[14]. We thus estimate the turbulent outer scale L of the 
solar wind stream to be 6.5 xlO 5 km, or ^1026 s in time 
units, compared with L = 0.35 in the MHD simulation. 

Increments 5B(r),5u(r) are considered to be “large” 
for our criterion when their magnitudes both exceed 1.5 
of their rms values. Using this threshold, we identify 
possible reconnection events in both datasets. See com¬ 
plete catalogues in [14]. Many of the candidate events 
in both datasets resemble the “double-step” magnetic 
reversals bounded by near-discontinuities, which Gosling 
tends to select with his original criteria [7]. However, 
we also see events with more gradual reversals over 
inertial-range scales in both datasets. In Fig. 1 we show 


2 



0.0 0.2 0.4 0.6 0.8 1.0 


0.0 0.2 0.4 0.6 0.8 1.0 


x/L 


x/L 


FIG. 1. Top: Event from the MHD simulation, at point 
(2.964,0.908,5.841) along a cut in the ^-direction. Magnetic 
fields (left) normalized by B rrns = 0.24, and velocity fields 
(right) by local upstream Alfven velocity Va = 0.7. Bottom : 
Event from Wind spacecraft data, on January 14, 2008, 13:50 
hr, normalized by B rrns = 2.5 nT, Va = 75 km/s. Distance 
x is normalized by L. MVF components are identified as L 
(red, solid), M( green, dotted), N (blue, dashed). 


events of this latter type. The vectors have been rotated 
into the minimum-variance-frame (MVF) of the magnetic 
field [ 5], calculated over the reversal region. The veloc¬ 
ities here (and in all following plots) are in a frame mov¬ 
ing with the local mean plasma velocity. Both events are 
inertial-range scale, occupying an interval of length 0.1 
in the MHD simulation and 2-3 min in the solar wind 
case. Although they do not have a “double-step” mag¬ 
netic structure, these two events do show the features 
characteristic of magnetic reconnection. There appears 
to be a reconnecting field component and an associated 
Alfvenic outflow jet in the L-direction of maximum vari¬ 
ance. A weak inflow is seen in the TV-direction of mini¬ 
mum variance, which is usually interpreted as across the 
reconnection “current sheet.” The M-direction of inter¬ 
mediate variance is nominally the guide-field direction, 
which in both events appears rather weak and variable. 

The MHD event shown in Fig. 1, top panels, arises 
from passage of the sampled ID cut close to a large, 
helical magnetic flux-rope appearing in the simulation. 
The maximum field strength in the rope is 8 times 
the rms strength in the database. Plotted in Fig. 3 is 
the original ID spatial cut, the magnetic cloud, and 
nominally incoming and outgoing field-lines along the bl¬ 
and L-directions of the MVF. There is a clear magnetic 
reversal, with incoming lines in the flux-rope twisting 
clockwise and into the page, but incoming lines to the 
left pointing out of the page. The field-line geometry 
is, however, quite complex since the lines exhibit the 
stochastic wandering assumed in the Lazarian-Vishniac 



FIG. 2. B-isosurface at half-maximum value 1.11 in 
B-lines sampled along N-direction in green and L-directions 
in red. The original ID spatial cut is the thick black line. 


theory of turbulent reconnection [ ]. Fig. 2 and all other 
spatial plots in this letter are available as 3D PDF’s [F ]. 

To identify large-scale geometry, it is necessary to spa¬ 
tially coarse-grain (low-pass filter) the magnetic field. 
For the theoretical basis of this coarse-graining approach 
to turbulent reconnection, see [ 6, 17]. Here we apply 
a box-filter with half-width L to obtain coarse-grained 
fields B, u, from which all inertial- and dissipation- 
range eddies are eliminated. The nature of the database 
event as large-scale reconnection becomes more evident in 
Fig. 3, which plots the lines of B. A central “X-point” at 
(2.84, 1.31, 5.73) was located by eye and a new MVF cal¬ 
culated in a sphere of radius L around that point. (This 
frame is rotated by ~ 20° in all three directions relative to 
the MVF for the original ID cut, but furthermore the M- 
and V-directions are exchanged). Field-lines are plotted 
at regular intervals along the L- and N -axes through the 
point. The plasma flow is incoming along the TV-direction 
and outgoing along the L-direction, and the magnetic 
structure is clearly X-type, with length ^ 0.4-0.6 (In¬ 
direction) and width ~ 0.15-0.2 (V-direction). Recon¬ 
nection events observed in the solar wind also appear to 
be X-type [7], although this structure has generally been 
interpreted in terms of Petschek reconnection. 

In fact, the field-line geometry of B in the MHD event 
is more complex than a single “X-point.” The complete 
structure is revealed by calculating the “perpendicular 
squashing factor” Q j_, a quantity devised to identify field¬ 
lines with rapidly changing connectivity in the solar pho¬ 
tosphere and corona [18]. We consider the Q_L-factor for 
the field-lines of B which begin and end on a sphere of 
radius 1.61 around the nominal “X-point” in Fig. 3. The 
Q ±-isosurface in Fig. 4 reveals a “quasi-separatrix layer” 
(QSL) whose cross-section has a clear V-type structure. 
The “hyperbolic flux-tube” (HFT) extending along the 
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FIG. 3. Same as Fig. 2, except for B rather than B. The iso¬ 
surface is for the half-maximum value |B| = 0.301. The blue 
vectors on the field-lines are u (in the local plasma frame). 

centers of these X’s has length about 0.54 and is aligned 
approximately with the M-direction, to within about 
35°. An HFT is the modern version of an “X-line” for 
3D reconnection, which does not usually admit true sep- 
aratrices and X-lines, and an HFT has the same obser¬ 
vational consequences as an X-line. It is thus interesting 
that very large-scale reconnection events in the solar wind 
(above integral scales) appear to have very extended X- 
lines, based on observations by multiple spacecraft [ L9] . 

Careful examination of the dynamics of this MHD 
event verifies that it is indeed magnetic reconnection, 
and fundamentally influenced by turbulence. The mag¬ 
netic flux-rope and the associated QSL persist over the 
entire time (0 to 2.56) of the database, drifting slowly 
with the plasma. The QSL and MVF also slowly rotate 
in time, with the MVF directions rotated through total 
angles ~ 40° at the final time and also the M- and N- 
directions exchanged around time 2.0. The time required 
for a plasma fluid element in the reconnection region to be 
carried out by the exhausts with velocities ~ 0.3-0.4 also 
happens to be about 2.0. Despite the high conductivity of 
the simulation, standard flux-freezing is violated in this 
event due to the turbulent phenomenon of “spontaneous 
stochasticity”, as we now verify. The exact stochastic 
flux-freezing theorem for resistive MHD [20] (which gen¬ 
eralizes ordinary flux-freezing), states that field-lines of 
the fine-grained magnetic field B are “frozen-in” to the 
stochastic trajectories solving the Langevin equation 

dx/dt = u(x, t) + V2X drf(t ), (1) 

where A = rjc 2 /An is magnetic diffusivity and rj(t) is a 3D 
Gaussian white-noise. The many “virtual” field-vectors 
B which arrive to the same final point must be averaged 
to obtain the physical magnetic field B at that point. We 
have chosen a point xj in the outflow jet in the +L-dir- 



FIG. 4. The quasi-separatrix layer Q± = 32 in cyan, and its 
cross-section in a plane normal to the M-direction in magenta. 


ection at time tf = 2 and solved (1) backward in time 
to to = 0, to find the positions of the initial points 
whose magnetic field vectors arrive at (xj,t/). This en¬ 
semble of points, plotted in Fig. 5, is widely dispersed in 
space. This disagrees with the predictions of standard 
flux-freezing, which implies that the ensemble should be 
close to a single point. In the lower panel of Fig. 5 is plot¬ 
ted the mean-square dispersion of this ensemble perpen¬ 
dicular to the L-direction, (rj_) as a function of reversed 
time r = tf —t. Consistent with previous results [6], the 
(backward) growth of perpendicular dispersion is diffu¬ 
sive (rj_(r)) ~ 8 At for very small r but then becomes 
super-ballistic, due to turbulent Richardson dispersion. 
As argued in [ ], the perpendicular spread in the time to 
exit with the outflow, y/ (rj_(2)) ~ 0.19, is close to the 
width of the reconnection region. This zone has both the 
width and the turbulent structure proposed in [4] , as can 
be seen also in Fig. 5 which plots in green the isosurfaces 
of the fine-grained current magnitude at half-maximum. 
There is a spatial distribution of many current sheets 
rather than a single large current sheet, as in laminar 
reconnection, and none of the sheets is located precisely 
at the QSL shown in Fig. 4. See [14] for 3D plots. 

The breakdown of standard flux-freezing is one evi¬ 
dence of reconnection in this event [21]. We have also 
verified that there is topology change of the lines of both 
fine-grained and coarse-grained magnetic fields. To show 
this, we decorate initial field-lines of either B or B at time 
t = 0 with a sequence of plasma fluid elements and then 
follow each element moving with the local velocity u for¬ 
ward in time to t = 2. We find that the plasma elements 
which initially resided on the same line at t = 0 end up 
on distinct lines at time t = 2 and some of these lines 
are outgoing in the +L-direction and others in the — in¬ 
direction. For movies, see [14]. We have also determined 
the average reconnecting electric field E rec for the large- 
scale magnetic field B, using a voltage measure proposed 
in [22]. We find that E rec ~ O.Olv aB in terms of local up¬ 
stream values va and B. Furthermore, at the length scale 
L of B, most of E rec is supplied by turbulence-induced 
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FIG. 5. Top: Vorticity isosurfaces at the half-maximum value 
uj — 69.1 in green. In red, origin points at time t — 0 of mag¬ 
netic field at final point (x/,t/) = (3.31,0.083,6.07,2.00). 
Bottom: Backward mean-square dispersion (t*j_(t)) orthogo¬ 
nal to L-direction as blue line. Reference curve (r*j_) = 8 At 
in green, and = 0.02r 8//3 in red. 

electric fields and resistivity gives only a tiny contribu¬ 
tion, always more than an order of magnitude smaller. 
These and many other detailed results for this event will 
be presented elsewhere. One finding is that this inertial- 
range event is not only highly 3D but also non-stationary 
in time. While the outflow jets are quite stable over time, 
the inflow is “gusty”, with variable magnitude and direc¬ 
tion veering in the N-M -plane (so that it is often in the 
nominal guide-field or M-direction). It is thus difficult to 
define an operationally meaningful “reconnection speed”. 

The main purpose of this Letter has been to present 
an example of inertial-range reconnection in MHD tur¬ 
bulence, to clarify its observational signatures. While a 
fluid description is surely applicable only to scales much 
larger than plasma micro-scales (e.g. the ion gyroradius 
in the collisionless solar wind), our simulation is remark¬ 
ably successful in reproducing observed features of large- 
scale solar-wind reconnection, together with crucial tur¬ 
bulent effects supporting theoretical predictions in [ , 5] . 
The characteristics are expected to change with length- 
scale, e.g. reconnection deeper within the inertial-range 
should have stronger guide-fields/smaller magnetic shear- 
angles [ , 23]. We also note some differences between the 
current MHD database and the solar wind, as our simu¬ 
lation is incompressible and isothermal, whereas the so¬ 
lar wind is slightly compressible and reconnection events 
there (including that in Fig. 1) often show enhancements 
of proton density and temperature in the reconnection 
zone. Furthermore, our MHD simulation has no mean 
magnetic field and is close to balance between Alfven 


waves propagating parallel and anti-parallel to field-lines, 
whereas the solar wind has a moderate mean field and the 
high-speed stream studied in this work is dominated by 
Alfven waves propagating outward from the sun. The in¬ 
fluence of these differences should be explored in future 
work seeking to explain large-scale solar wind reconnec¬ 
tion in detail within an MHD turbulence framework. 
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